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In the Letter the initial conditions (IC's) of cosmological N-body simulations by the 
Virgo Consortium [5] are analyzed and it is concluded that the density fluctuations are rather 
different from the desired ones. We have repeated the analysis of the IC's using our own code 
and the code provided by the authors of 1 , obtaining results that disprove the criticisms. 
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Fig. 1 - Left: density variance in balls of radius r (in units of the box size), estimated with our code 
and the one provided by the authors of Right: estimated power spectrum. The resolution limit 
corresponds to a distance « 2(A). 



Fig'-GJshows the density variance, <T^{r), computed for the initial conditions of the SCDMl 
model of the Virgo Consortium |23 . For (r) we used the estimator quoted in T : variance 
of density fluctuations in 20000 spheres centered at random positions. The theoretical CDM 
behaviors have been computed with the power spectrum employed by the Virgo Consortium 
12], and they correct the wrong ones in [T]. For completeness, we also show the estimated 
power spectrum, P{k). We find that the estimated value of cr^(r) does follow the theoretical 
CDM prediction in an intermediate range of scales 2(A) r < L/2 ((A) being the mean 
interparticle separation and L the box size). On scales outside this trustworthy range, the 
expected systematic departures due to discreteness and finite-size set in. In the same way, 

© EDP Sciences 



EUROPHYSICS LETTERS 



0.004- 0.006 



resolution limit 



analytical _ 

VIRGO ICs 




i 




— 





i 0.010 0.012 



resolution limit 



analytical 



. own code 

code of [ 1 ] 




r/box 



Fig. 2 - Absolute value of the estimated correlation function; the inset panel shows the oscillatory 
behaviour of ^(r) for small separations. Left: random balls are allowed to overlap. Right: they are 
constrained to avoid overlap. 



the estimate of P{k) also agrees with the desired CDM behavior in the intermediate range 
27r/i < fc < 7r/(A). 

Fig. 121 shows the estimated £,{r). We used the estimator quoted in in turn a simple 
modification of the one for cr^ (variance in spherical shells centered at randomly chosen par- 
ticles). The estimate of ^(r) exhibits noticeable fluctuations around zero, whose amplitude 
is larger than the desired CDM decay. This amplitude is comparable to the variation of the 
estimate when using different random seeds for the random sphere centers. Such a noise also 
arises when estimating ^(r) for simple test cases whose correlation is known theoretically (un- 
correlated particles, particles at the nodes of a lattice). This suggests that these fluctuations 
are mainly noise of the ^-estimator employed in In these circumstances, a systematic de- 
parture from the CDM behavior provoked by the method how the IC's are generated cannot 
be assessed. 

We also learnt that the results shown in were derived in reality with a modification 
of the quoted estimators: the probing spheres were constrained to avoid mutual overlap, so 
that their number is effectively smaller than the intended 20000 for r large enough. We have 
observed that this leads to no detectable significant difference when estimating a, but the 
estimated ^(r) tends to be larger in the large-r end compared to the unconstrained estimate. 
This sensitivity on the number of probing spheres throws further doubts on the reliability of 
the measured £,(r), and it suggests looking for better estimators. 

Our conclusions are: (i) cr'^{r) and P{k) reproduce very well the desired CDM behavior 
in the expected intermediate range of scales; (ii) the estimator of ^(r) employed in ^ is too 
noisy to allow one to claim that ^(r) does or does not follow the expected behavior; (iii) hence, 
these measurements provide no evidence against the procedure employed to generate the IC's. 
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